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SOME EFFICIENT METHODS FOR OBTAINING INFINITE SERIES SOLUTIONS 
OF n** 1 -ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS 

by Gabriel Allen 
Lewis Research Center 

SUMMARY 

ir. 

In this report, some methods of obtaining series solutions for n - order linear 
ordinary differential equations are presented. The required analysis of the differential 
equation to determine whether the point of expansion is an ordinary point or a regular 
singular point is shown. The standard method of obtaining series solutions is presented 
next and, for regular singular points, the procedures are catalogued according to the 
value of s, the difference between the two solutions to the indicial equation. Theorems 
concerning the number of series solutions to second-order linear ordinary differential 
equations are given. Examples are worked out for each possibility. 

Next, the definition and some properties of the theta operator are shown, and this 
is followed by a description of the use of this operator in obtaining the indicial equation 
and recurrence relation for n -order differential equations. The method is then ap- 
plied to the same examples that were worked out by the standard method in order to 
show the reduced effort involved in using the theta operator. 

Finally, the generalized hypergeometric function is defined and some relevant prop- 
erties described. These properties are then applied to the same examples described 
previously. 

A short reference section is included which lists all the needed properties of gen- 
eralized factorials, theta operators, and generalized hypergeometric functions for ob- 
taining series solutions to differential equations. 


INTRODUCTION 


Most practical problems in science and engineering are described in terms of re- 
lations between two or more physical quantities. Such relations are often in the form 



of n th - order linear ordinary differential equations. Solutions to such equations can 
rarely be expressed in closed form. Usually, the solutions must be obtain numerically. 
Even in those cases where closed-form solutions can be obtained, numerical methods 
must almost always be used in order to obtain specific M y versus x" curves or tables. 

Therefore, the method of solution by power series has great utility, since the form 
of the solution is ideally suited for obtaining a "y versus x" table even when the closed 
form of the solution is not known. Furthermore, series solutions are ideal for com- 
puting machines, which would otherwise often resort to them anyway to compute tran- 
scendental functions given to them in "closed form” solutions. 

The preceding list of advantages of series solutions notwithstanding, a fair amount 
of mathematical manipulation is usually necessary in order to obtain an explicit series 
solution (i.e. , to obtain easily computed expressions for the coefficients of the general 
term in the series). 

This report points out that there exist some methods for greatly reducing the amount 
of labor usually involved in such tasks. Sadly enough, these methods are not new, but 
seem to have been neglected in this application even by otherwise sophisticated and 
well-informed applied mathematicians. 

In this report, some preliminary definitions and theorems concerning the series 

J.V 

solutions of n -order linear ordinary differential equations. of complex variables are 
listed. The theorems cover the existence of solutions about ordinary points and regular 
singular points. Following this section, the "standard" method of Frobenius is re- 
viewed. A good treatment of this method is shown in reference 1. Next, the use of the 
theta operator is demonstrated. Finally, the use of generalized hypergeometric func- 
tions is described. Some illustrative examples thread through the entire treatment and 
are solved by each of the methods. The superiority of the theta- operator and 
generalized-hypergeometric-function methods over the "standard” method is made 
patently clear, thereby. 

The following two points, which are also discussed in the relevant part of the text, 
should be noted: 

(1) Some of the definitions and all the theorems pertain to second-order linear 
ordinary differential equations only. 

(2) All the examples shown have two- term recurrence relations. The theta- 

i-U 

operator method, however, can be applied to n - order and will also yield recurrence 
relations with an arbitrary number of terms. 

It should be emphasized that neither the theta operator nor the generalized hyper- 
geometric functions are, in any sense, an invention of the author. His contribution is 
in the emphasis and application of these two mathematical entities towards reducing the 
pain associated with obtaining series solutions to linear ordinary differential equations. 
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Readers already familiar with the method of Frobenius and the analysis of singular 
points according to the solutions to the indicial equation can go directly to the section 
THETA-OPERATOR METHOD (p. 20). 

SYMBOLS 

n-1 

(a) generalized factorial symbol, TT (a + i) 

n . i=0 

pFq(z) generalized hypergeometric function with p numerator parameters and 

q denominator parameters 

D d/dz, differentiation operator 

R(c) value of real part of complex number c 

s difference between two roots of indicial equation of a second- order linear 

differential equation; taken so that R(s) s 0 

9 z(d/dz), theta operator 


SOME DEFINITIONS AND THEOREMS FOR SECOND-ORDER EQUATIONS 


The material in this section applies only to second-order differential equations. 


Definitions 

Standard form . - The standard form for a linear second-order ordinary homoge- 
neous differential equation is defined to be 

+ p(z) — + q(z)w = 0 (1) 

dz 2 dz 

Ordinary point . - The point Zq (in the finite plane) is an ordinary point of equa- 
tion (1) if and only if both p(z) and q(z) are analytic at z Q . 

Singular point. - If Zq is not an ordinary point of equation (1), it is a singular 
point. 
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Regular singular point . - Let P(z) = (z - Zq)p(z), and Q(z) = (z, - q(z). Then 

the point Zq is a regular singular point of equation (1) if and only if 

(1) Zq is not an ordinary point. 

(2) Both P(z) and Q(z) are analytic at Zq. 

Irregular singular point . - If Zq is a singular point of the differential equation (1) 
but is not a regular singular point, it is an irregular singular point. 

Point at infinity . - To determine the nature of the point at infinity, equation (1) is 
transformed by using the new independent variable £ = 1/z and examining the resulting 
equation at £■= 0. The transformed differential equation is 


? 4 A, + 2{ 3 dw _ { 2 p /l\dw + /l\ w = „ 

d |2 d? U/d« W 


In standard form, this becomes 


+ 7 r(|) — + /c(|)w = 0 (la) 

as 2 d « 


where 




77 (|) = 


2|- P 



and 



Therefore, a description of the point z = °° is determined by the analysis of the be- 
havior of 77 ( 1 ) and k(!) at £ = 0. 

As an example, examine the following hypergeometric differential equation for sin- 
gular points: 


z(l - z) + (2 - 3z) — - w = 0 


dz 


dz 
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In standard form, this equation is 

d 2 w | 2 - 3z dw _ 1 w _ 0 

dz 2 z(l - z) dz z(l - z) 

Here, p(z) is not analytic at z = 0 or z = 1 so that both points are singular points. 

o 

Here zp(z) = (2 - 3z)/(l - z) is analytic at z = 0, and z q(z) = -z/(l - z) is also an- 
alytic at z = 0. Therefore, z = 0 is a regular singular point. Similarly, 

(1 . z) p( z ) = !^?z 
z 


is analytic at z = 1 and 


(1 - z) 2 q(z)= - LLl-Zl 
z 

is analytic at z = 1. Therefore, z = 1 is also a regular singular point. 

To examine the point at °°, the transformation £ = 1/z is made, the differential 
equation (1) is rewritten in the standard form (la), and the point £ = 0 is examined. In 
this case, 

1 

4 

and 


1 

« 2 

so that £ = 0 (and, hence, z = °°) is a regular singular point. 

An example of an irregular singular point is the point at °° of the differential equa- 
tion 


^ + k 2 w = 0. 


dz 


k = a constant 
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The transformed equation becomes 


+ 2 | 3 dw +k 2 w = 0 


dr 


d£ 


which in standard form is 


d 2 w 



2 dw 

I d£ ^4 


w = 


0 


Clearly, x(|) fails the test for a regular singular point, and thus z = °° is an ir- 
regular singular point of the differential equation. 


Theorems 

The following two theorems are useful in determining the existence of series solu- 
tions for second- order linear ordinary differential equations. 

Theorem 1. - If Zq is an ordinary point of the differential equation (1), then there 
exists a series 


w 


•S’- 

n=0 


n( z ' z o) n 


( 2 ) 


which satisfies the differential equation within the circle of convergence of the series, 
has a nonzero radius of convergence, and has two arbitrary constants a Q and a^. It 
should be recalled that the radius of convergence R is at least the distance from the 
point of expansion z^ to the nearest singularity. 

Theorem 2. - If Zq is a regular singular point of the differential equation (1), then 
there exists a series 


Wj(z) 


S s(’ - 


n=0 


, n+c. 


(3) 


which satisfies the differential equation within the region of convergence of the series 
and has a nonzero radius of convergence. In equation (3), c^ is defined as follows: Let 
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Cj and c 2 be the roots of the indicia! equation (see eq. (7)), and let R(c^) and R(c 2 ) 
be their respective real parts. Then c^ is chosen so that RCc^) > R(c 2 ). 

STANDARD METHOD OF FROBENIUS 

While a knowledge of the usual method is not needed for a discussion of the two non- 
standard methods described in the next two sections, there are, nevertheless, two use- 
ful purposes to be served by reviewing this method before introducing the others. First, 
since it is already familiar to the reader, it is a pedagogically efficient way to lay the 
framework for series solutions to differential equations. Secondly, a proper apprecia- 
tion of the difference in algebraic labor involved can only be obtained by a direct com- 
parison of several methods of solution to the same differential equation. In the examples 
to be shown, therefore, the first method of solution is the method of Frobenius. The 
other methods are then compared with it. 


Solution About an Ordinary Point 

The main reason for including even a brief discussion of these solutions in this re- 
port is to emphasize the difference between solutions about an ordinary point and solu- 
tions about regular singular points. The following properties of solutions about ordinary 
points should be emphasized: 

(1) There always exist two linearly independent series solutions about every ordi- 
nary point of a second- order linear ordinary differential equation. 

(2) There is no indicial equation for solutions about an ordinary point. The solution 
is a power series in integral powers of z - Zq. 

By comparison, in the case of solutions about a regular singular point, theorem 2 
only guarantees that one of the two linearly independent solutions will be in the form of 
a power series. The second solution may also be in this form, but sometimes it cannot 
be expressed in a power series. Secondly, the power series solution is usually in frac- 
tional, rather than integral, powers of z - Zq. If c^ = 0, the series solution in the 
case of regular singular points will also be in integral powers of z - z^; but this is the 
exception, rather than the rule. These distinctions are referred to again in the ex- 
amples. 
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. Solution About a Regular Singular Point - The Method of Frobenius 

It is convenient for this discussion to write equation (1) in the form 

P(z) dw + QW „ = 0 (lb) 

dz 2 z - z 0 dz (z - z Q ) 2 

If Zq is a regular singular point of the differential equation (1), both P(z) and Q(z) are 
analytic at Zq and we may write 

OO 

P(z) = ^ P n ( z " z o) n i z " z 0 i < R i ( 4 ‘) 

n^O 

Q(Z) = £ %( z ~ z oT l z " Z 0 I < R 2 (5) 

n=0 

where Rj and R2 are the radii of convergence of the two series. 

The first part of the procedure consists of assuming a solution of the form 

w = E a n( z - z o)" +C (6) 

n=0 

where a Q ^ 0 and is otherwise arbitrary and c and the other a n 's are to be deter- 
mined. The series (4), (5), and (6) are then substituted into equation (lb). The result 
is 

J (n + c)(n + c - l)a n (z - Zg) n+C 2 + J P m ( z " z o) m Y (n + c)a n( z " z 0 T +C 2 
n=0 _m=0 n=0 

+ fe q i( Z - Z o)fE a n( z - z or C - 2 ]=° 

j=° JL n=0 

The lowest power of z - Zq in this equation is c - 2. Equating the coefficient of this 
power to zero results in the equation 
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c(c - l)a 0 + P Q ca 0 + q Q a 0 = 0 


By hypothesis, ag / 0. Invoking this hypothesis results in the important indicial equa- 
tion of the differential equation (lb) 

c 2 + c(p 0 - 1) + q Q = 0 (7) 

In order to continue the discussion, it is essential to define the quantity 

s = c t - c 2 (8) 

where c^ and Cg are the two solutions to the indicial equation (7) and R(c^) s R^). 
There are three possibilities, each governing a different procedure. 

Case a: s / an integer or zero . - In this case, two power series solutions exist, 
one for each solution to the indicial equation. The two solutions will therefore be of the 
form 


w 0 


w 


1 


oo 

z 

Tl = f 


I 


n=0 



(9) 


where ag and bg are arbitrary. 

Case b: s = 0 (or c^ = Cg). - In this case only one solution to equation (lb) exists in 

power series form. The second linearly independent solution is of logarithmic form. 
Therefore, the two solutions are 


w n = 


I a n(* - 

n=0 


,n+c, 


J 0) 



> 


(a logarithmic form) 


( 10 ) 
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Case c: s = a positive integer . - The procedure to be followed in this case is to 
/ \ n+c 2 

begin in powers of (z - Zq\ , where Cg is that solution to the indicial equation with 
the smaller real part. 

The trial series is substituted into equation (lb) and, beginning with the lowest 
power, the coefficients of each power of z - Zq are combined and then set equal to 
zero. In this way, some of the a n 's are determined. This procedure is continued until 
the coefficient a g is reached. At this stage of the procedure, it is determined if the 
equation has been a "lucky" example or an "unlucky" one. 

A lucky case is one in which a is found to be arbitrary, and two series solutions 

s 

exist. On the other hand, in an "unlucky" situation, it is found that an attempt to deter- 
mine a leads to a contradiction, that the differential equation has only one series 
solution, and, furthermore, that series solution which does exist is the one for the other 
solution to the indicial equation, so that all the work to this point has been wasted. The 
second solution in the unlucky case is again a logarithmic solution. 

Therefore, with the aid of theorems 1 and 2 and an examination of the roots of the 
indicial equation, it is possible to determine in advance, except for case c, how many 
power series solutions exist and to estimate the difficulty involved in trying to obtain 
these solutions explicitly. 

The following examples will serve to clarify much of the procedure described pre- 
viously. 


Examples 

Example 1: Solution about an ordinary point . - Consider the series solution about 
z = 0 for the differential equation 


w" + 2zw’ - 8w = 0 


(ID 


This differential equation is already in standard form and z = 0 is an ordinary point of 
the given differential equation. Therefore,, by theorem 1., two. solutions of . the form 

oo CO 

L a Z n will exist and both a n and a., will be arbitrary. When / , a z is substi- 
n 0 1 n= 0 n n 

tuted into the differential equation, the result is 


Z n(n_ 1 } v n ' 2 

n=0 


+ £ 2na n z n - £ 8a n z" = 0 
n=0 n=0 
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or 


^ n(n - l)a n z n " 2 + ^ (2n - 8)a n z n = 0 
n=0 n=0 

- 2 - 1 

The two lowest powers of z are z and z , which result in the requirements 

0 x a Q = 0 

0 x a^ = 0 

so that both a^ and 0 a^ are arbitrary, in accordance with theorem 1. Shifting index 
n — n + 2 in the first sum, combining the two resulting sums, and equating the coeffi- 
cient of z n to zero result in the recurrence relation 

(n + l)(n + 2)a n+2 = (8 - 2n)a n n > 0 


or 

n(n - l)a n = (12 - 2n)a n _ 2 n a 2 


When aQ is used as the arbitrary constant, only even powers result. When n = 2k, 
the recurrence relation becomes 


a 2k = 


12 - 4k 

2k(2k - 1) 2k_2 


2(3 - k) 

k(2k - 1) 2k ' 2 


k > 1 


Therefore, 


= a .= 4a 


1 ( 2 - 1 ) 


0 ” 0 


_ 2(3 - 2) „ _ d 2 _ 4 Q 
a a 2 

4 2(4 - 1) 1 3 3 U 


a 2k ~ 0 k ~ 3 
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so that Wq(z) = Sq|i + 4z 2 + (4/3)z^j. 

When is taken as the arbitrary constant, only odd powers (n = 2k + 1) result. 
Therefore, in this case, 


_ 5 - 2k k>l 

215+1 "k(2k + l) 2k - 1 

This series does not terminate, and the general expression in terms of a^ is found to 
be 


a _ (5 - 2)(5 - 4). . , (5- 2k) 

215+1 k! [l x 3 x 5 . . . (2k + 1)] 1 

o 

Therefore, the second solution (which is linearly independent of Wq(z)) is 

(5 - 2)(5 - 4) ... (5 - 2k) z 2k+1 
1 x 3 x 5 ... (2k + 1) k! 

Example 2: Solution about a regular singular point where s / 0 and s / an in- 
teger. - Consider the series solutions about z = 0 of the differential equation 

2zw” + w' - w = 0 (12) 

It is clear that z = 0 is a regular singular point of this differential equation. Writing 
the differential equation in the form (lb) we obtain 


w^z) = 


ajZ 


+ a. 


Z 

k=l 


w” 



z 

w = 0 


(12a) 


The first requirement in a series solution about a regular singular point is the indi- 
cial equation. Here, pg = 1/2, qg = 0 so that equation (7) becomes 



Therefore, c = 0 and 1/2 and s = 1/2. This means that the differential equation be- 


12 



is sub- 


longs to case a and there will be two series solutions, 
stituted into equation (12), the result is 


When w = Y] a z n+c 
n=0 n 


^ 2(n + c)(n + c - l)a n z n+c-1 + ^ (n + c)a n z n+c_1 
n=0 n=0 


oo 

-E 


a n z 


n+c 


n=0 


= 0 


The indicial equation has already been obtained, so it is not necessary to equate the 
coefficients of the two lowest powers to zero individually. Therefore, a direct attack 
on the recurrence relation can be made at this point. For this purpose, the index is 
shifted on the first two sums and the resulting equation becomes 

[(n + c + l)(2n + 2c + l)a n+1 - a n Jz n+c = 0 

n=0 


The recurrence relation is 


n+1 


(n + c + l)(2n + 2c + 1) 


n > 0 


or 


n 


ct 1 

n-1 

(n + c)(2n + 2c - 1) 


n > 1 


(12b) 


This recurrence relation can be used for both solutions to the indicial equation. How- 
ever, the symbol for the coefficients should be changed in the second solution to avoid 
confusion. 

The allowed values of c are 0 and 1/2. They are substituted one at a time into 
equation (12b) to obtain the two linearly independent solutions to equation (12). When 
c = 0 is used, 


a n-l 

n(2n - 1) 


and the general term becomes 
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Therefore, 


a 0 2 n 
^ — — 

n n'.[l x 3 x 5 . . . (2n - 1)] (2n)l 0 


w Q (z) = a Q / z n 

° °/ ^2nV. 

n=0 


For the second solution (c = 1/2), it will be less confusing if the symbol for the coeffi- 
cient of z n+ ^ 1 / 2 ^ is denoted by b n instead of a n> Then equation (12b) becomes 


n 


n+ l\(2n) n(2n + 1) 

v 2/ 


The general term in this case becomes 


Therefore, 


2V 


b = 

n (2n + 1)! 


oo 



2 n z n +(V2) 

(2n + I)! 


There are a number of clues which distinguish w^(z) corresponding to c = 0 from 
a solution about an ordinary point. For differential equations which have two-term re- 
currence relations, the solutions about an ordinary point always break up into a series 
of even powers and another one of odd powers, there being two solutions of integral 
powers. In this case, Wq(z) has all powers of z, both even and odd. Furthermore, of 
course, the second solution is a series of fractional powers of z, rather than integral 
powers as would be the case for a solution about an ordinary point. 

Example 3: Solution about a regular singular point where s - 0 . - If a slight change 
is made in the preceding example, the case s = 0 results. The differential equation is 


zw" + w' - zw = 0 


(13) 
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Clearly, z = 0 is a regular singular point of this differential equation; Pq = 1 and 
qQ = 0. Therefore, equation (7) is 


c 


2 


- 0 


In this case, c = 0 and 0 and s = 0. Only one solution about z = 0 exists in power 
series form for this differential equation. The procedure for obtaining this solution is 

OO 

E n+c 

a n z into equa- 
tion (13) and combining terms of like powers of z, the result is 


n+c+1 


Z / \2 n+C -1 \ A 11-fL, t*J. ft 

(n + c) a n z - 2a V ■ 0 
n=0 n=0 


The index is now shifted in the first sum and the following recurrence relation is ob- 
tained: 


n 


ri+2 


(n + c + 2 )‘ 


n > 0 


or 


n 


a n- 2 
(n + c) 2 


n > 2 


The coefficient a< has been required to be equal to zero by virtue of equating the co- 
0 1 

efficient of z to zero, so that no odd values of n will appear in the solution. The 
allowed values of n can therefore be most readily obtained by setting n = 2k, and the 
general term becomes 


k 

JJ ( 2 i + c ) 

j=i 

Therefore, the series solution to equation (13) is 


15 






,2k+c 


T 2 I 


(13a) 


T[ (2j + c) 
3=1 


j 


Now 


■ w ° = ( w 'L 


or 


w 0 (z) = a Q 


1 + 


QO 

I 


,2k 


— 


7i s2k ( k! ) 2 


Xtf 

(k :) 2 


k=0 


The second solution is 


w l = 



The procedure is standard but somewhat lengthy. Since the methods to be shown in 
the subsequent sections would not be helpful in reducing the effort involved in finding the 
second solution, the procedure is not exhibited in this report. Nevertheless, the solution 
itself is shown so that the description of the solution as not being expressible in series 
form will be fully appreciated. Thus, 


oo 



/ 

k=l 


(13b) 


where 
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is the sum of the first k terms of the harmonic series. Both terms on the right side of 
equation (13b) diverge as z — 0. However, the sum of the two terms remains finite. 

Example 4: Solution about a regular singular point where s is a positive integer - 
lucky case. - Consider the differential equation 

9z 2 w M + 3z(z 2 + 2)w' + (z 2 - 2)w = 0 (14) 

In standard form, the equation becomes 


,, z 2 + 2 , z 2 - 2 n 

w" + w' + w = 0 


3z 


9z 


Therefore, z = 0 is a regular singular point with Pg = 2/3 and qg = -2/9. Under these 
conditions, equation (7) is 



Thus, c = 2/3 and -1/3 and s = 1. 

The procedure when s is an integer is to try using the smaller of the solutions to 
the indicial equation, that is, Cg. The reason for trying Cg is that it is then found that 
either the recurrence relation for the coefficient a g is impossible or else it will be 
automatically satisfied, leaving ag and a g arbitrary and yielding the general solution 
from the work involving Cg alone. When the recurrence relation for a g leads to a set 
of inconsistent conditions, no solution using Cg is found. The inconsistent condition 
which occurs is the requirement that ag be 0 - contrary to the hypothesis that ag be 
the coefficient of the first nonvanishing term of the series and, therefore, must be arbi- 
trary. In such a case, the only series solution to the differential equation is one in- 
volving Cp rather than Cg. In this example, two solutions are found by using Cg. 

Substitute w = V] a into equation (14): 

JrTn n 


I[4 -3HH - -3- • I •(■ ■ 

n=0 n=0 
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or 


£ 9n(n - l)a n z n -< 1 / 3 > + £ 3na n z n+ < 5 / 3 > = 0 
n=0 n=0 

By equating the coefficients of z -1 / 2 and z 2//2 to zero, the following requirements 
result: 


0 x ap = 0 
0 x a 1 = 0 

so that ag and a^ are both arbitrary. Note that s = 1, so a g = a^ here. The recur- 
rence relation is found to be 


a = - n -~-— ■ a 9 n > 2 

n 3n(n - 1) n ' 2 


Therefore, in this case, a^ =0, k>l, and 


w 0 ( z ) = a Q z 


-1/3 


(14a) 


The solution involving a^ is obtained by letting n = 2k + 1 and using 


2k - 1 

l 2k+ 1 6k(2k + 1) a2k_1 


The general term is then found to be 


L 2k+1 


(-if 


6 K k! (2k + 1) 


a^ k > 1 


The second solution, which is also a series solution, is then 
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w 1 (z) = aj 


QO 

I 


(_l)k z 2k+(2/3) 
6 k (2k + l)k! 


It should be noted that the theorem guaranteeing a series solution about a regular sin- 
gular point specifically refers to using c^ in the solution. In the "lucky" case, the 
second series solution always involves c^ since it is a solution whose lowest power of 
z is Cg + s, which is c^. 

Example 5: Solution about a regular singular point where s is an integer - unlucky 
case . - A common differential equation which illustrates this case is Bessel’s equation 
of index 2: 

z^w" + zw’ + (z 2 - 4)w = 0 (15) 

The usual analysis shows that z = 0 is a regular singular point and that Pq = 1 and 
qg = 4, so that the indicial equation is 

c 2 - 4 = 0 

Here c = ±2 and s = 4. As in the previous example, c 9 = -2 is tried first. There- 

oo " 

E n _ 2 

a z is substituted into equation (15). After simplifying, the results 
n=0 n 
can be written as 

2 »<" - 4 >v n ' 2 + 2 v n = ° 

n=0 n=0 


Since s = 4, we will equate coefficients of successively higher powers of z to zero until 

_ o 

we have come to a^. The coefficient of z is 

0 x a Q = 0 

so aQ is arbitrary. The coefficient of z“* is -3a^ = 0, so that a^ = 0. The coeffi- 
cient of z° requires contributions from both sums now: 

2(-2)a 2 + ag = 0 


so a 2 = a 0 /4. Next, 
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-3a 3 + = 0 


o 

so a 3 = 0. Finally, we come to the coefficient of z : 

0 x a^ + a 3 = 0 

This requires that ag = 0. But a 3 = Sq/ 4 from a previous requirement. Therefore, 
the only way to keep the complete set of requirements on the coefficients consistent is 
to set ag = 0 also. This requirement is the one which is inconsistent with the hypothe- 
sis that a.Q be arbitrary. The result is that no series solution using c 3 = -2 can be 
found . 

n i O 

The only series solution for Bessel's equation is a series in powers of z . The 

second solution is found by a method somewhat similar to the one for the s = 0 case, 

but the algebra is often still more involved. The solution is expressed as w 3 (z) = 

(3w /3c\ and is also a logarithmic solution. Neither solution will be displayed ex- 

\ ' c=c 2 
plicitly in this report. 


THETA-OPERATOR METHOD 

Two procedures which are used in obtaining series solutions for differential equa- 
tions are often lengthy. These are (1) finding the indical equation in solutions about regu- 
lar singular points and (2) finding the recurrence relation. Admittedly, equation (7) does 
not require much effort compared to the usual "textbook” method of equating to zero the 
coefficients of the lowest power of z - Zq. However, the theta operator yields the in- 
dicial equation and the recurrence relation automatically. As will be seen from the ex- 
amples, the theta operator greatly reduces the total effort needed to obtain series solu- 
tions . 

Furthermore, the method is applicable to n 1 order, not just to second-order 
linear ordinary differential equations. This method has not been shown in any widely 
used textbooks, although it has been known for many years. Recently, however, Brand 
has finally included a good section on it in his book on differential and difference equa- 
tions (ref. 2). Since he does not include the "standard" method in his book, the reader 
cannot appreciate the method as much as the comparison shown herein will permit. 
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Definitions and Properties 

The definition of the theta operator is 

6 = z — =zD (16) 

dz | 

where D = d/dz. 

Since we are interested in n -order linear ordinary differential equations we must 
be able to express D n in terms of 6. To begin with, we show that 

e(e-i) = z 2 b 2 ( 17 ) 

Proof : 

0(0- 1) = z— (z — - l\ - z /z— +— V — = zV Q.E.D. 
dz \ dz / \ dz; 2 dz J dz 

Next, the generalized factorial notation is introduced: 

n-1 

(a) = a(a + 1). . . (a + n - 1) = fT (a + j) (18) 

n j =0 

In terms of the T function, (a) n = r(a+n)/T(a). By using equation (18), the following 
relation can be proved by induction: 

z n D n = (-l) n (-0) n (19) 

Proof : The relation certainly holds for n = 1 (by definition). The second part of the 
proof, therefore, consists in showing that, if it holds for n = j, it must hold for 
n = j +1, j > 1. Now, if equation (19) holds for n = j, then 

zW = (-l) j (-0). (19a) 

Operate on equation (19a) by (-1)(- 0 + j). The right side becomes (-1)^ +1 TT (- 6 + j), 

i+1 k=0 

which equals (-1) J (-0)j + ^. Next, the left side is examined: 
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e(z^) •= zD(z^) = z(jz^ - 1 D^ + Z y +1 ) = jz^D^ + z^ +1 D^ +1 
Since (-1)(- 6 + j) = (8 - j), the left side becomes 

(0-j)(z^D^) = jz^D^ + z^ + ^D^ + ^ - jz^D^ ■= z^ + ^D^ + ^ 

Therefore, z^ + *D^ + * = (-l)^ + '*'(- 6). + ^, j >1, and the second part of the proof is complete. 

The superiority of the 0-operator to the D-operator arises from the effect of the 
former on z 11 : 0z n = nz 11 . Therefore, any polynomial expression in 8 operating on 
z n yields the same polynomial in n times z 11 . This property is useful enough to be ex- 
pressed as a numbered equation. Thus, 


P(6)z n = P(n)z n (20) 

In equation (20), the left side indicates a differential operation on z n , whereas the 

right side is an ordinary algebraic expression. 

An additional property of the 8 - operator should be mentioned. It is sometimes con- 

k 

venient to make a change of independent variable of the form u = az . After making such 
a change of variable, the operator u d/du is more convenient to use than 9 itself. If 
this new operator is denoted by <p, then 



( 21 ) 


Probf: 


u 


az 


du 


k dz d 
du dz 


az 


k 1 


d _ k 

— cLZ 


du dz 
dz 


kaz 


J- = l z A. 

k-1 dz k dz 


8 

k 


Note that for the differential equation in <p-form, the new s is 1/k times the old one. 


Application of Theta Operator to Series Solutions of n* h -0rder 
Linear Ordinary Differential Equations 


A 


general n 


th 


-order linear ordinary differential equation is usually written in the 


form 
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|c 0 (z)D n + c 1 (z)D n ' 1 + . . . + c n (z)Jw = 0 


We shall limit ourselves to differential equations for which all the coefficients c^(z) are, 
at most, polynomials in z. Then, by using equation (19) whenever it is needed, such a : 
differential equation can always be put into 6-form. The result is 

[a o ( 0) + zA^fl) + . . . + z m A m (6)]w = 0 (22) 

where all the A. (6) are polynomial operators in 6. ■ 

In order to show how the 6-operator can be used to solve such a general differential 
equation, we assume that z = 0 is a regular singular point so that w is of the form 

OO 

X) a z n+c . Upon substituting this series for w in equation (22) and using property (20), 
n=0 n 

there results the equation: 


J a n A 0 (n + c)z n+c + £ a^n + c)z n+c+1 + ...+£ a n A m (n + c ) zll+C+m = 0 
n=0 n=0 n=0 

(23) 

Clearly, the smallest exponent of z is c; and, in order that a Q be arbitrary, the fol- 
lowing indicial equation is automatically obtained: 

A 0 (c) = 0 (7a) 

The recurrence relation will be obtained once the index has been shifted by an ap- 
propriate amount in each series after the first one. Upon shifting the indices so that the 
sum in each series is over z n+c , equation (23) becomes 


°o 

Yj [ A 0 (n + C)a n + A l (n - 1 + c)a n- 1 + • ’ ' + A m (n " m + c)a n-m] (z " Z 0 )n+C = 0 
n=m 

(24) 

The recurrence relation (which is an n-term one) is therefore 

A Q (n + c)a n + A x (n - 1 + c)a n _ 1 + . . . + A m (n - m + c)a n _ m - 0 (25) 
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The preceding results were obtained for a solution about z = 0. If a solution about some 
Zq ^ 0 is desired, the following simple modifications will make the procedure for 
Zq / 0 equivalent to that for Zq = 0. 

A transformation to a new independent variable z’ = z - Zq is made. Then each 
polynomial coefficient c^z) becomes a polynomial in z’, that is, c?(z').. Finally, a 
new operator 6 ' = z'(d/dz') is defined. After all these changes are made, an equation 
whose form exactly parallels that of equation (22) results, with Aj(6’) and z’ playing 
the same roles in the new equation as A^fi) and z play in equation (22). 

It is, therefore, clear that there is no loss of generality in using an expansion about 

z 0 = 9 * 

The following steps summarize the 6-operator method for obtaining series solu- 
tions. to n th -order linear ordinary differential equations: 

(1) Write the differential equation so that equation (19) can be applied. 

(2) Write the differential equation in 6-form (i.e. , eq. (22)). 

(3) If a solution about a regular singular point is desired, the indicial equation is 
A q (c) = 0. 

(4) Write the recurrence relation using equation (25). 

Of course, the rest of the problem, such as obtaining an explicit formula for a n , must 
be done as previously. Naturally, the analysis which determines the nature of the point 
about which the solution is to be taken must precede the application of the 6- operator 
method. The following examples allow a comparison to be made of the amount of effort 
necessary to obtain solutions by the 6-operator method and the ’’standard” method. 


Examples 

2 2 

Example 1: 2z(d w/dz ) + (dw/dz) - w = 0 . - It was shown in the previous section 
that z = 0 is a regular singular point, so the method can be applied immediately. In 
order to write this differential equation in 6-form (see eq. (19)), the equation must first 
be multiplied by z. The resulting equation is 


2z 2 + Z - zw 


dz 


dz 


0 
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Now the conversion to 6-form becomes trivial: 

[2 0( 0 - 1 ) + 0 - z]w = 0 


or 


[6(26 - 1) - z]w = 0 

Here, A Q (0) = 6(26 - 1), and A^6) = -1. Therefore, the indicial equation is c(2c- 1) = 0, 
which, of course, has the same solutions as were obtained by the standard method of 
Frobenius. 

The recurrence relation is obtained immediately. It is 

(n + c)[2(n + c) - l]a n - a n-1 = 0 


or 


n (n + c)(2n + 2c - 1) 

This is the same as equation (12b), and the procedure is unchanged from the "standard” 
method from here on. 

Example 2: zw" + w* - zw = 0 . - As in example 1, this differential equation must 
be multiplied by z before equation (19) can be applied. After this has been done, the 
result is 


2 2 
z w" + zw'-z w = 0 


or, in 6-form, 


|<9(0 - 1) + 6 - z 2 ]w = 0 


which simplifies to 


(e 2 - z 2 )w = o 
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2 2 
Therefore, A g(0) - 0 , A^(0) = 0, and A 2 ( 0 ) = -1. The indicial equation is e = 0, and 

the recurrence relation is (n + c) 2 a„ - a„ 0 = 0. Since the only series solution for this 

11 11” & / O 

case (s = 0) is for c = 0, the relation a n = a n _ 2 / n > where n > 2, is obtained as be- 
fore. 

2 2 2 

Example 3: 9z w" + 3(z + 2)zw’ + (z - 2)w - 0 . - This equation is ready to be 
converted to 0-form as it stands. The result is 

[90(0 - 1) + 3(z 2 + 2)0 + (z 2 - 2)]w = 0 

which simplifies to 

[90 2 - 30 - 2 + z 2 (30 + l)]w = 0 

Here, Aq(0) = 90 2 - 30 - 2 = (36 - 2)(30 + 1), A^(0) = 0, and A g(0) =30+1. The indicial 
equation is (3c - 2)(3c + 1) = 0, which results in s = 1, as before. The recurrence re- 
lation is [3(n + c) - 2][3(n + c) + i]a n + [3(n + c - 2) + l]a n _2 = 0, n > 2. The rule in this 
case is to try the c = -1/3 solution. This attempt leads to the recurrence relation 


a 


n 


3n - 1 - 6 + 1 a 
(3n - 1 - 2)(3n -1 + 1) n “ 2 


3n(n - 1) n ~ 2 


n > 2 


which was obtained by the standard method of Frobenius. 

Of course, in working out a case where s is an integer, the coefficients of powers 
of z must be separately equated to zero until the coefficient a g has been shown to be 
arbitrary. There is no point in looking for a recurrence relation until after the exist- 
ence of a series solution for c = Cg has been established. The 0- operator method can- 
not eliminate this part of the problem, but putting the differential equation into the form 
of equation (23) makes it easier to carry it out. In this example, equation (23) would 
take on the special form 


J 9n(n - Da,/-' 1 / 31 * £ 3na n zMV3) + 2;L 0 


Therefore, the coefficients of the two lowest powers are readily obtained from the equa- 
tions: 0 x ag = 0, and 0 x = 0. Once it has been shown that these two coefficients 
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are arbitrary, the recurrence relation can be found with confidence that the solution for 
c = Cg exists. 

It is hoped that these three examples have demonstrated the. advantages of using the 
^-operator in obtaining series solutions to linear ordinary differential equations. 


SOLUTIONS USING GENERALIZED HYPERGEOMETRIC FUNCTIONS 

It turns out that, for most of the linear ordinary differential equations which have 

n+c 

found extensive application, the coefficients of the general term z can be written 
without resorting to a recurrence relation at all. Furthermore, all the series solutions 
that exist for these differential equations can be written with their coefficients com- 
pletely determined. This is so because all these functions satisfy differential equations 
which are special cases of a single differential equation (27). 

Admittedly, there do exist some linear ordinary differential equations which are not 
of the form permitting this to be done, and, for these, the use of the theta operator is 
recommended to reduce the labor involved in obtaining series solutions. However, a 
great many differential equations in common use do lend themselves well to the method 
described in this section, and series solutions for these equations can be obtained with 
a minimum of effort. 


Definitions and Properties of p Fg(z) 

Definition. - The definition of a generalized hypergeometric function as used herein 
is (see eq. (18)) 


oo 



The a^ and b^ are called numerator parameters and denominator parameters, re- 
spectively. It should be noted that the n! in the denominator is not included as a pa- 
rameter in the notation. Nevertheless, as will be shown later, its presence cannot be 
ignored in constructing additional series solutions. 
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Another notation in common use is more cumbersome than the left side of equa- 
tion (26) but gives more explicit information. This alternate notation is 

F(a 1 ,. . .,a p ; b p . . .,b q ; z) 

In the alternate notation the subscripts p and q are optional. The "grammar" in this 
notation is important. The list of denominator parameters is set off by semi-colons to 
denote the separation of these parameters from the numerator parameters to the left 
and from the variable to the right. 

When written out, the right side of equation (26) takes the form 

&i x a n x ... x a a« (a* + l)a«(ar> + 1) . . . a (a + 1) ~2 - 

1 + £ p z + -±-h £_£ P_p L. + . . v + 

bj x b 2 . . . x b q b l( b l + 1 ) b 2^2 + 1) • • • b q( b q + 1)2! 

From this form it is clear that if any a i = 0, p F q (z) = !• 

Convergence of F (z). - There are p distinct factorial factors in the numerator 

r 4 

and q + 1 such factors of n in the denominator (the n! must be included in consider- 
ations involving convergence). Therefore, applying the ratio test to equation (26) shows 
that pF q (z) converges everywhere if p < q + 1, converges only at z = 0 if p > q + 1, 
and converges inside the unit circle if p = q + 1. The ordinary hypergeometric function 
of Gauss is a 2 F 1 (z), and so it converges in the region |z| < 1. 

A few common examples of functions expressible as F (z) are 

r 4 



ln(l + z) = 


y tif^y ^. z y m ts£.. 

Z— * n + 1 Z— 4 (n + l)n! Z— t (2) n! 
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sin -1 z - z 2 F 1 (l/2,l/2; 3/2; z 2 ) 
tan _1 z - z 2 F 1 (l/2, 1/2; 3/2; -z 2 ) 


= 


n 


r(n + 1) 


F, (- ; n + 1; -z74) 


0 A 1 


Application of Generalized Hypergeometric Functions to 
Solutions of Linear n th -Order Differential Equations 


The complete set of generalized hypergeometric functions is too general for the 
applications. Since it is desirable to have convergence over a region more extensive 

v 

than the point z = 0, it is necessary that p < q + 1. Furthermore, if any denominator 
parameter were equal to a negative integer, -p, then the (p + l) th term in equation (26) 
would have a zero in the denominator. As shown by Rule I (p. 30), the same result will 
occur in a related p F^ if any two denominator parameters differ by an integer. 

Therefore, we will limit our discussion to generalized hypergeometric functions 
which are subject to these restrictions. Rainville used the ©-operator to derive the fol- 
lowing differential equation which has such generalized hypergeometric functions as 
solutions (ref. 3): 


e ]”f (© + b. + 1) - z JJ (0 + a. 
3=1 J i=l 


w = 0 


m 


where b^ is not a negative integer and no two b's differ by an integer. 

This differential equation is very convenient to use since the denominator param- 
eters are unambiguously separated from the numerator parameters. The convenience 
of this form is due entirely to the fact that the differential equation is in ©-form; such 
convenience disappears when the differential equation is written in a conventional manner 
(see example 1, p. 32). 

By applying to equation (27) the discussion in the two preceding main sections, the 
following facts can be written at once: 

(1) The point z = 0 is a regular singular point. 
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(2) The indicial equation is 


q 

c Jf (c + b. - 1) = 0 
j=l J 


Therefore, the q + 1 solutions to the indicial equation are c = 0 (common to all p F q' s ) 
and c = 1 - bj, where j = 1, . . . , q. 

(3) A 0 (e) = 6 fr (6 + b - 1); A 1 (e) = - J[ (6 + a.) 
u jL'l 3 1 i=l 1 

Therefore, the recurrence relation is a two-term one. Since, in this section, the sym- 

bol a^ is reserved for the i Ln numerator parameter, let a n be the coefficient of z n 
in the series for Wg(z). Then the recurrence relation takes the form 


a = 
n 


Jl^ (n + c - 1 + a t ) 


a , n > 1 
n-1 


(n + c) J[ (n + c - 1 + b.) 
i = l J 


st 

The linear ordinary differential equation (27) is of (q + 1) order and, therefore, 
has q + 1 linearly independent solutions. One of these is 


w 0 (z) = p F q( a p- • *’ a p ; b l»‘ • • ’ b q’ z ) ( 28 ) 

The other q linearly independent solutions can be obtained from Wg(z) by using the fol- 
lowing rule: 

Rule I: Add 1 - b m to every parameter in the numerator and the denominator (including 
the (l) n disguised as n!), and multiply the resulting function by z ^ m . 

Thus, 


1-b 


w m ^ z ) = z p F q^ a l + 1 ' b m’ • • ’ ’ a p + 1 " b m ; b l + 1 " b m’ * • ' ’ b m-l 


+ 1 - b m- 1 + 1 - b m’ b m+l + 1 ‘ b m V 1 ~ V z > < 28a > 


th 

The n! thus becomes the rn denominator parameter, 2 - b . The new _F_ is of 

1 ’ m p q 
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the proper form, since it also has an n! in the denominator obtained from tHe param- 
eter which was in w n (z). This parameter becomes b_ + 1 - b_ = 1 in the new 
m 0 mm 

pV 

Remarks Concerning Reduction of Second-Order Linear Ordinary 
Differential Equations to the Form of Equation (27) 

It is recommended that the theta operator be used in all instances, so that it is as- 
sumed that the differential equation is already in the form of equation (22). The fol- 
lowing remarks may then be helpful in reducing second-order linear ordinary differential 
equations to the form of equation (27): 

Remark (1) . - Equation (27) indicates that a differential equation whose solution is a 
pFq has only a two-term recurrence relation. Therefore, unless the differential equa- 
tion in form (22) has a two-term recurrence relation, it cannot be put into the form (27), 
and no solution of the form p Fq exists. 

Remark (2) . - If equation (22) has a two-term recurrence relation, it will be of the 

form [a o ( 0) + z k A k (0)Jw = 0. At this point, a change of independent variable t = -z k 

should be made. This transformation will result in a new equation: [A 0 (k<p) - 

tA.(k<p)]w = 0, where cp = t d/dt. Following this transformation, the equation should be 
K 2 

divided by the coefficient of cp in A p and the new polynomial should be factored, re- 
sulting in (cp - c^)(cp - c 2 ), where c^ and c 2 are the solutions to the new indicial 
equation in which t is the independent variable. At this stage, a decision must be made 
as to whether a change of dependent variable will be useful. 

Remark (3) . - In making this decision, it may be noted that the theta operator has the 
following effect on a function w = z n v: 

6(z n v) = z n (n + 0)v (29) 

This property may be used to transform any term Aq(6)w of the form (6 - c^)(0 - c 2 )w 

Cj C| c< 

to z a A 0 (e + c^)v, which is of the form z 1 e{e + (cj - c 2 )]v, where w = z A v. It 
should be noted that s still equals c^ - c 2 in the new differential equation having v as 
the dependent variable. Therefore, the following useful rule has been established: 

C 1 

Rule II : Given that A p (0)w = (8 - c^)(0 - c 2 )w, then the transformation w = z v will 
result in a differential equation in which this term will be replaced by 6{8 + s)v. 
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Rule II enables a simple analysis to be made of solutions in the s = 0 case. 

P q 2 

Wheh s = 0, the rule shows that Aq(6)w will be replaced by 6 v. Therefore, if a 

pFq is a solution to this differential equation, the (single) denominator parameter will 

have to be b = 1. Therefore, b - 1 = 0, and the application of Rule I shows that there 

will never be a ’’second" F solution about a point for which s = 0. 

P q 

If s is a positive integer, q, then Aq(6)w will transform to 8(6 + ju)v, and b will 
equal 1 + ju. Therefore, if a F solution exists, Rule I should apply for finding the 

' ir H. 

second solution. 

Remark (4) . - The discussion in remark (3) has been limited to an examination of 
Aq( 0). It is certainly a necessary condition that Aq( 8) be of the proper form in order 
that F be a solution to the differential equation. However, it is not a sufficient con- 
dition. It is also necessary for A^( 8) to be of the proper form. 

See examples 4 and 5 in the next section to appreciate how these remarks can be ap- 
plied. 

Examples 

Example 1: The hyper geo metric equation. - In order to gain some confidence in the 
method, it will first be applied to an example which is tailormade for it, the hypergeo- 
metric equation of Gauss: 

z(l - z)w" + Jy- (a + (3 + l)zjw’ - apw = 0 

The equation must be multiplied by z in order to get it into the 6-form. The resulting 
equation is 

(1 - z)z^w" + ^y - (a + 1 3 + l)zjzw’ - a |3zw = 0 
In 6-form, the equation becomes 

j(l - z)6(6 - 1) + Jy - (a + + l)zje - Q?j3z|w = 0 

or in the form of equation (27) 

[e( 6 + y - 1) - z(6 + a)(6 + /3)jw = 0 
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Therefore, one solution to this differential equation is 


w 0 (z) 


2 f 1 («,/3; r; z) 



(a) (B) n 
v V H/ n z 

(y) n n: 


The second solution is obtained by adding 1 - y to each parameter, including the 1 

1- v 

which represents n! , and multiplying the new 2 F i by z • The result is 

Wj(z) = z 1_y 2 F i( a + 1 - y,0 + 1 - y; 2 - y; z) 


Thus, the two linearly independent solutions can be written almost on sight with all the 
coefficients of z n in the series completely determined. 

Now that we have seen how to use the method in a tailormade case, some of the ex- 
amples used in previous sections of this report will be examined next. 

Example 2: 2zw" +w' - w = 0 . - This differential equation in 6-form was shown 
to be 


je(20 - 1) - zjw = 0 


In order to put this equation into the form of equation (27), the equation must be divided 
by 2. This will allow the parameters to be correctly determined. Finally, then the 


form 


X® ■ i) - 


tion which is a so 


|w = 0 results. In this form, it is clear that the hypergeometric func- 
lution to the differential equation has one denominator parameter and 
no numerator parameters and that the independent variable must be z/2, rather than z. 
Therefore, 


w 0 ( z ) = 0 F 1 (-; 1/2; z/2) 

Furthermore, the second solution is 

w x (z) = z 1_(1/2) q F 1 (~ J 1 + 1 - 1/2; z/2) = z 1/2 gF^-; 3/2; z/2) 

The solutions in this form have completely determined coefficients and are certainly 
adequate for any computations which may be desired. In our case, however, we may 
feel more comfortable if we show explicitly that these are, indeed, the same solutions 
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that were obtained previously. This task is readily accomplished with the use of the 
following relation which is listed in the next main section: 



For instance, in series form, 



If this is substituted into the series for Wq(z), the result is 


w 0 (z) 


Z 2V z^ 

(2n)I if. 


which agrees exactly with what was obtained by the standard method of Frobenius. 
Similarly, 
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Now, 




3 

2 


+ n - 


1 

2 


(2n + l/A) = (2n ± 1): 

W„ 2 2 V. 


Therefore, 

1 _ 2 n 
2 n n :/- N j (2n + 1)! 

\ 2 / n 

and Wj(z) also agrees with the result obtained by the standard method. 

Example 3: the s = 0 case (case b, p. 9) . - The example of case b, for which 
s = 0 can be used to illustrate a second- order linear ordinary differential equation which 
has only one as a solution. In 6-form, the equation is 

( 8 2 - z^)w = 0 

Since this is an s = 0 case, only one series solution will exist. As recommended in 

o 

remark (3), the change of independent variable t = z should be made, thus transform- 

2 2 

ing the differential equation to (4 cp - t)w = 0, or, upon dividing by 4, (<p - t/4)w = 0. 

This equation is in the form of equation (27), and the first solution is a qFj with 
t/4 as the variable: 


CO 



n=0 


This is the same solution as was obtained by the standard method of Frobenius. How- 
ever, the rule for obtaining the second solution fails here, since b-l=l-l = 0, so 
the ’’second" solution is the same as the first. Even hypergeometric functions cannot 
give a second series solution where none exists! 
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Example 4: the lucky s = integer ease (case c, p. 10). - This differential equation 
in conventional notation was 9z2w M + 3(z^ + 2)zw’ + (z2 - 2)w = 0. In the section 
THETA-OPERATOR METHOD, it was shown that it could be written in 6-form as fol- 
lows: 


[(36 - 2)(36 + 1) + z 2 (36 + l)]w = 0 
2 

The change of independent variable t = -z is indicated and results in the transformed 

differential equation: (6 <p - 2)(6 cp + 1) - t(6© + 1) w = 0. Upon dividing by the coefficient 
2 

of cp , this equation becomes 



An application of rule II is indicated. Therefore, the dependent variable is changed by 
using w = t 1//3 v to result in the differential equation 



This equation has a solution of the form 

v Q (t) = 1 K 1 (l/2; 3/2; t/6) 


from which one may obtain 

w 0 (z) = t 1/3 v 0 (t) = (-z 2 ) 1/3 1 F 1 (l/2; 3/2; -z 2 /6) = (-l) 1/3 z 2/3 ^(1/2; 3/2; -z 2 /6) 

This is the solution which was obtained as w^(z) by the standard method of Frobenius, 
To reconcile the forms, we may use equation (26) 


z 2 / 3 


1 F 1 (l/2; 3/2; -z 2 /6) = z 
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Using the property (a + 1) = (a +n)/a (a) n given in the compendium (p. 40), then, 
if a = 1/2, J 



Therefore, (l/2) n /(3/2) n = l/(2n + 1). Also, (-z 2 /6) = (-l) n z 2n /6 n . These changes 
convert the solution obtained here to the exact form of the solution obtained by the stand 
ard method of Frobenius. 

If Rule I is applied to find the second solution to the differential equation here, the 
result is 


Vl (t) = f l/2 1 F 1 (0; 1/2; -z 2 /6) = f 1/2 x 1 
according to the note following the written-out series from equation (26). Therefore, 
W 1 (z) = tV3 t - 1/2 = t -l/6 = ^. z 2y 1/6 = ( . ir l/6 z -l/3 

which agrees with Wq(z) in the section STANDARD METHOD OF FROBENIUS. 

Example 5: the unlucky s = integer case - Bessel’s equation . - Consider Bessel’s 
equation of order n, where n is a positive integer: 

z 2 w” + zw' + (z 2 - n 2 )w = 0 
In 0-form, this equation becomes 


(e 2 - n 2 + z 2 )w = 0 

This is the case where s = 2n and is a positive integer. Following the recommended pro- 

2 

cedure in remark (3), the change of independent variable t = -z converts this equation 
to 


(4 <p 2 - n 2 - t)w = 0 
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which, upon dividing by 4, becomes 




For this equation, c^ = n/2, Cg = -n/2 and s = n. Next, the change of dependent vari- 
able, w = t n / 2 v, is made. The differential equation with v as the dependent variable is 



This equation is of the form (27) with the denominator parameter b = n + 1. There- 
fore, Vg(t) = o F i( - ; n + 1; t/4). From this solution, 

w Q (z) - t n / 2 v = z n Q F 1 (-; n + 1; -z 2 /4) = 2 n T(n + l)J n (z) 


from the short list of common functions expressible as p F q' s - Thus, the series for 
J n (z) has been obtained in a relatively few lines. As we saw before, the second solution 
to Bessel’s equation is a logarithmic one. Since the denominator parameter is n + 1, 

1 minus this parameter equals -n, and the rule for obtaining the second solution would 
result in a qF^ with a denominator parameter of 1 - n. Since denominator parameters 
are not allowed to be zero or negative integers, a second solution in terms of hypergeo- 
metric functions cannot be obtained for integral n. It should also be noted that if n = 0, 
then b = 1 and the second solution is the same as the first (this is an s = 0 case). 

It is hoped that these examples are sufficient to show the use of the theta operator 
and hypergeometric functions in obtaining series solutions to n -order linear differen- 
tial equations. 


SUMMARY OF PERTINENT APPROACHES TO OBTAINING SERIES SOLUTIONS 
OF LINEAR ORDINARY DIFFERENTIAL EQUATIONS 

The first step in obtaining a series solution to a linear ordinary differential equa- 
tion is an analysis of the point about which the series is to be expanded (see section 
SOME DEFINITIONS AND THEOREMS FOR SECOND-ORDER EQUATIONS). Such an 
analysis gives information concerning the number and kinds of series solutions which 
exist about that point. Following this analysis, it is recommended that the differential 
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equation be put into theta form. At that point an attempt to put the equation in the form 
(27) may be made. If this attempt is successful (and the bj’s satisfy the restrictions 
preceding equation (27)), all q + 1 solutions may be written with the coefficients com- 
pletely determined. 

In those relatively rare cases in which the solutions cannot be (readily) expressed 
as pFq' s > the indicial equation and recurrence relation can be written for whatever so- 
lutions are expressible in series form. 

The following section is a compendium of the formulae and properties of theta opera- 
tors and generalized hypergeometric functions. It should be complete enough to enable 
these procedures to be followed without reference to external sources. There will prob- 
ably be exceptions to this, but the author believes them to be rare. 
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COMPENDIUM OF PROPERTIES OF THETA OPERATOR AND 
GENERALIZED HYPERGEOMETRIC FUNCTIONS 


Generalized Factorial Notation 


n- 1 


(a) n = TT (a + j ) = F ^ a + - -- (see eq. (18)) 
j=0 r(a) 


(a) n+1 = (a + n)(a) n 


( a + l) =l^ = ^hn ( a) 
a a 


(a) 


2n 


2 / n 2 2n/a ± l) 


For the special case where a = 1, 


( 1 ) 


2n _ (2n): _ 1x3x5. . . (2n - 1) 

r,n 


2/ n 2 2n (l) n 2 2n n! 


This formula yields a simple expression for the product of odd integers - a fairly com- 
mon occurrence. 


Theta Operator 


The theta operator is defined 


6 = z ~ (see eq. (16)) 
dz 


z n D n = (-l) n (-e) n (see eq. (19)) 
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If P(0) is a polynomial in 6, then 


P(0)z n = P(n)z n (see eq. (20)) 

An n - order linear ordinary differential equation in theta form is 

Ja q ( e) + zA 1 (6>) + . . . + z m A m (e)]w = 0 (see eq. (22)) 
The indicial equation is 

Aq(c) = 0 (see eq. (7a)) 


The recurrence relation is 

A Q (n + c)a n + Aj(n - 1 + c)a n _ 1 + . . . + A m (n - m + c)a n _ m = 0 (see eq. (25)) 

k 

For a change of independent variable - if t = az , then 9 = k cp (see eq. (21)), 
where cp - t d/dt. In the differential equation in cp- form, the solutions to the indicial 
equation are c 1 /k and c 2 /k. 

For a change of dependent variable, 

0(z n v) = z n (0 + n)v (see eq. (29)) 


Therefore, 


(6 - c 1 )(e - c 2 )w - 6(6 + s)v (see Rule II) 


where w = z v. 


Generalized Hypergeometric Functions 


The definition and notation is 


pFq( z ) F ( a p- • -’ a p5 bp- • -jkqJ z) - 



n=0 L J 


3 

T ( a i) 

i=l n , 




n: 


(see eq. (26)) 


The differential equation for written in theta form is 


P 


1 TT (0 + b. - 1) - z TT (0 + a.) 

j=l ] i=l 1 


w = 0 (see eq. (27)) 


The q + 1 solutions to this equation are 


w 0 (z) = p F q (a l’* • *»V b l»- * -’V Z) 


W (z) = z 
m 


1-b 


m 


V a l +1 " V- • ” a p + 1 - b m ; b x +1 - b m ,. . .,b m _ 1 


1 " b m’ 2 " b m’ b m+l + 1 " b m’ • 1 * ’ b q 


+ 1 - b m ; z) m = 1, . . . , q (see eq. (28a)) 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, August 7, 1972, 

503-10. 
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